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Abstract 

We investigate the CGHS model through numerical calculation. The behavior 
of the mass function, which we introduced in our previous work as a "local mass" , 
is examined. We found that the mass function takes negative values, which means 
that the amount of Hawking radiation becomes greater than the initial mass of the 
black hole as in the case of the RST model. 
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1. Introduction 



Two dimensional dilaton black hole models have attracted a great deal of atten- 
tion in recent years [P], [2], [3], f|, [5], [| [7], They provide a way to treat black holes 
that evaporate due to quantum effects and so may serve as a useful toy model for 
real four dimensional black holes, whose evaporation processes lead to such issues 
as the well-known information problem. Although the analysis of 2d dilaton black 
hole models has so far remained within the semi-classical approximation, this is 
still a great step towards the complete understanding of the quantum behavior of 
black holes. 

We are interested in studying the end point of these black hole models within 
the semi-classical approximation. Numerical studies of the CGHS model have been 
made by several authors || [T0|, 11|, but, little is known about the end point of 2d 



black holes. In particular, it is a difficult task to estimate the total amount of 
Hawking radiation produced by a black hole in a numerical study. 

In this paper we investigate 2d dilaton black hole models with the local mass 



introduced in our previous paper [JL^]. This local mass can be related to the integral 
of the energy momentum tensor in the weak field region, and hence we can use it 
to calculate the total energy radiated from the black hole. 

In the next section, we review our previous work and introduce the local mass. 
We will use this method to investigate the analytical solution for evaporating black 
holes in section |3]. In particular, we demonstrate that the excess of Hawking 
radiation over the initial incoming energy leads to a negative value of the local 
mass. 

We present the result of our numerical study in section [TJ It is shown that the 
CGHS model suffers from the same problem, i.e., the excess Hawking radiation, as 
in the RST model. Section [)] is devoted to a discussion. 



2. Local mass 



In this section we will explain the local mass which we have used in [[12] to 
analyze 2d quantum black hole models. An expression of the local mass for spher- 
ically symmetric four dimensional gravity can be found in the literature (see e.g. 
|13||). Tomimatsu [H| used the local mass in four dimensions to investigate an 
evaporating Schwartzschild black hole. The local mass for the 2d dilaton black 
hole is also considered in Jl5|, [16], |17| . 



The gravitational part of the action which we are going to study consists of the 
2d metric g, a dilaton and a constant A: 



2tt J 



d 2 x^g e- 2<t> (R + A(V<pf + 4A' 



(2.1) 
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For the matter part, we introduce N scalar fields ff. 



12.2) 



Throughout the paper we restrict our investigation to initial matter configurations 
such that a shock wave of fi comes in at an advanced time Xq with the total energy 
M. 

In the Hamiltonian formalism, we see that a certain combination of the two 
Hamiltonian constraints for the gravitational part gives the total divergence of the 
local mass [E|. Decomposing the two dimensional metric g as 
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(2.3) 



one can rewrite the gravitational action in terms of the canonical variables with 
the Hamiltonian constraints: 



S (G) 



7T^0 + 7T 7 7 - N H - N{Ki . 
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It then turns out that the following combination of the above constraints becomes 
exactly the total derivative of a scalar function: 
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where the prime, stands for spatial derivative. The above scalar function M. can 
be expressed as 

-20 + Ae -20 



M 



(2.6) 



A A7 
We call the function M. the local mass. 

It would be more transparent to execute the procedure described above in the 



covariant formalism |[17|| . The energy-momentum tensor of the gravitational action 
is 



T (G) = 



2tt 5S^ 



-9 8g 
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If we make a current J„ from the energy-momentum tensor, 



(2.8) 
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it turns out that it is divergence-less, V M J7^ = identically. This is the reminiscent 



of what Kodama [18fl found in the four dimensional spherically symmetric system. 
Because of its divergence-less nature there exists a scalar function such that = 
g /j,v£ ua V a-M. ■ In other words, the following identity holds: 



The explicit form of M. is given by 



M 



-2(j> 



A 



(2.9) 



(2.10) 



which is the covariant expression of 

There is, however, an ambiguity in the above construction, as M. may be mul- 
tiplied by a constant or have a constant added to it. This ambiguity is fixed in the 
following way. We can relate the local mass to the flux of the energy momentum 
tensor in the weak field region. The flat space-time solution for this model is known 
as the Linear Dilaton Vacuum (LDV): 



ds 2 



-da + da . 



(2.11) 



It is easy to see that in this geometry ( 2.11|) the local mass Q2.10T ) becomes zero. 
For the above geometry the identity (|2.9|) simply becomes 



3 + m = -(r|? + t| g) ; 
d_M = t {g 2 + t {G }. 



(G) i 

Now T_j__ is nothing but the gravitational part of the equations of motion which 
should equal zero if quantum effects are negligible. Hence, for the flat space-time 
we obtain 



d_M = T 



r(G) 
++ > 



(GO 



(2.12) 



The ambiguity which may exist in the definition of the local mass is fixed by the 
above relations ( |2.12| ) and the fact that the local mass becomes zero for LDV ( |2.11[ ). 

If we take an asymptotically flat coordinate system, the relations ( |2.12|) also 
hold at infinity, where the gravitational field is weak. Assuming that the local mass 



-4- 



becomes zero at (x + = 0,x = — oo), integration of the first relation in ( |2.12 ) along 
past null infinity (x~ = — oo) starting from x + = gives 



M(x + ,-oo) = f (-T. 
Jo 

In particular, the local mass at (x + = +» 



do- + . (2.13) 

= — OO 

-oo) gives the ADM mass: 
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If we consider the integration of the second relation in ( p. 12] ) along the x direction 
at future null infinity we obtain 



M{+oo,x ) = M{+oo,-oo) + I T. 

- f (-t! g) ) 

J — oo 



(G) 



M 



x+=+oo 



da' 



(2.15) 



which coincides with the expression of the Bondi mass. 

The shock wave of fi, which represents the implosion of the matter shell, gives 
the following energy momentum tensor: 



T 



(M) / + - A 
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where Uq = e Xx o . From ( |2.13| ) we obtain 

M(a + ,a 



x oo 



M0(o + - <r+) . 
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The above expression agrees with the result from the direct calculation \V2\ of the 
local mass for the CGHS solution. 

We emphasize that the local mass is the amount of flux of the energy-momentum 
tensor in the weak field region. We will see that it reflects quantum effects. We 
can calculate the local mass on the matter shock-wave line 



as 



M\ 
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The local mass Ai diverges at the quantum singularity, where e~ 2</> reaches ^. 



(2.19) 
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It is worth noting that the local mass ( 2.10|) can be rewritten as 



A< = i-( c -a* - g) J2 , (2.20) 

where R stands for the scalar curvature of two-dimensional geometry. To derive the 
above expression (|2.20 ) we have used the equations of motion (4.1) which include 



the effective term of the quantum correction for the matter part. Then one can 



compare our local mass with another mass in the literature [ 10 1, 



M e// = ^(l-^)le- 2 ^. (2.21) 
Both of them coincide at future null infinity with 

m(aT) = lim -U" 2 ^, (2.22) 

x+— >oo 4A 

which is appeared in j|. 

We have introduced a function, "local mass", which should be useful to investi- 
gate 2d dilaton black holes. Unfortunately, incorporating quantum effects prevents 
us from obtaining an analytical solution for the system. Therefore, we shall study 
it numerically. On the other hand Russo, Susskind and Thorlacius Q found a 
model which is exactly solvable even with the quantum correction. In the next 
section we analyze the behavior of the local mass for the RST model. 



3. Behavior for RST model 

In this section we estimate the value of the local mass for the RST model. In 
Ref. Q it was demonstrated that incorporating an (artificial) effective term 



N 
48^ 



J d 2 Xy/=g<pR, (3.1) 



into the CGHS model, as well as a term from the quantum effect of the matter, 
one obtains the exact solution: 

, 24 _ 2 , 24A 2 a; + x- 1. ,. 2 , , 24M(x + - a^K , , + . 

* = N G N 2 HX X X ] M^ 6{X ~ Xt) ' 

(3.2) 

ds 2 = -e^dx+dx-. 

The above solution describes the formation of a black hole and its evaporation 
through Hawking radiation. 
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RSTmassa.eps 



Figure 1: The behavior of the local map "^r the RST model is shown in 3Dplot 




Figure 2: The behavior of the local mass for the RST model in density contour 
plot, where the density decreases as the local mass grows. 



We can visualize the RST model by the behavior of the local mass at each space 
time point. It is shown in Figures |I| and |2|. 

In the region x~ — ► — oo, the local mass shows the same behavior as that of 
the original CGHS model. It is clear from the fact that the above RST solutions 
reduce to the classical ones of the CGHS model at x~ — > — oo. On the other hand, 
the local mass behaves differently at 

M\ x+ _ x+ = M ~f X ° X ~ /AQ - (3.3) 
x ~ x o -\ 2 x£x- - N/48 

This difference is caused by the additional effective term ( |3.1| ). 

The RST model has the advantage that we can evaluate the value of the local 
mass at x + = +00 analytically. The information at future null infinity contains 
that of the end point of the black hole. Also, at future null infinity the local mass 
directly corresponds to the total energy radiated so far. The result at future null 
infinity is 

M\^(x-) = M + ^{-3^ + Alog(l - -3^)}. (3.4) 

The fact that in the RST model the black hole falls into negative energy state 
can be reproduced exactly in our study with the local mass. In Ref . [|J it is claimed 
that up to the point where the singularity and the apparent horizon meet together, 

, _ N\Xq , asm /NX _ 1 \ 

s - ASM [ h 

x ~ — — 1 r 35 ) 

the total amount of the energy radiated by the Hawking radiation becomes 

M+— (l-e-^M/NXy 

48 V ; 

This exceeds the energy supplied by the matter shock wave at x + = Xq by ^(1 — 
e -A&M/NXy j n f ac ^ l oca l mass at the future null infinity gives 

M(+oo, X :) = -^(l-e-*™n. (3.6) 

This confirms the relation between the local mass and the energy momentum tensor 
( p,15| ). In Figure one can see the region where the local mass becomes negative 
at the upper right of (xf,xj). 
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4. Numerical study of the CGHS model 



The equations of motion of the CGHS model at the one-loop level are given by 



2(1-^ e 2 ^) 9+«9_0 (4(9+0 9_0 + A 2 e 2 ?) = 



(4.1) 



9+«9_0 = (l - - e 2 ^ «9+<9_p , 

where p stands for the so-called Liouville field: 

ds 2 = —e 2p dx + dx~ . 

These equations which incorporate the quantum effects are no longer exactly solv- 
able and we shall solve them numerically. 

The parameter A can be scaled away and we put A = 1. The initial data 
surface, the line of a shock-wave of an infalling matter and we will 

choose Xq = 1. The and p on x + = 1 are those of the linear dilaton vacuum 
(LDV) and those values are given by 

0(1, x~) = p(l, x~) = - 1 - log (-x-) . (4.2) 

We shall impose boundary conditions at past null infinity such that the solutions 
become the classical ones there. Actually at a large negative value of x~, we impose 
that the and p take the classical values. 



For fixed x + the equations (|4.1| ) reduce to the ordinary differential equations for 
(9+0 and d + p with respect to x~ . We numerically integrate the equations from a 
large negative value of x~ through the apparent horizon, (9+0 = 0, to the singularity 
e 2(/> _ \2/ N. In practice we calculate up to the singularity. With the calculated 
values of 9+0 and <9+p we step forward in the x + direction. Actually the values 
of 9+0 and <9+p on the shock- wave line x + = 1 can be calculated from Eqs. (|4. 1|) 
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9-1-0(1, X 



d + p{l,x~ 



1 M 1 

2 2 ^Z^^J_ X - -N/12 ' 

(4.3) 




and they can be used to check the accuracy of the numerical algorithm. The 
numerical method which we adopted is explained in the appendix. 
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Figure 3: The singularity and the apparent horizon in the CGSH model. 



The result of the singularity and the apparent horizon in the CGHS model for 
M = 25 and N = 600 is shown in Fig. || It indicates that the singularity and the 
apparent horizon meet at finite x + and x~ which we denote (x% P ,x^ P ). 

The local mass function ( |2.6|) in the CGHS model for the same M and N is 
shown in Fig. [| and [5|. The figure |5| shows that the local mass becomes 

negative in some region, before the x direction and after the endpoint 

x + > x\ P + (some const.) in the x + direction, in a manner similar to the RST 
model. 



5. Discussion 



In the previous section we found that the amount of the Hawking radiation in 
the CGHS model exceeds the incoming energy as was seen in the RST model case. 
In other words, the CGHS model fails to end up as the vacuum but evolves into a 
space-time whose scalar curvature is negative. 

One interpretation is that this merely means the failure of the CGHS model 
as a toy model for four dimensional black hole models. Note in particular that we 
have coupled scalar matter fields to the gravitational field in a way that is different 
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CGHSmassa.eps 



Figure 4: The local mass function in the CGSH model (3Dplot). 




Figure 5: The local mass function in the CGSH model (density plot). 



from the four dimensional case. This coupling of the matter fields enables us to 
calculate the back reaction but could also cause the excess of the radiation energy. 

One possible modification of the CGHS model may come from string theory. 
There have been many attempts to resolve the problems associated with black 
holes by applying string theory and one approach along this line is the following. 
The term {l/4)e- 2 ^a'R abcd R abcd ^, where a' is the inverse of the string tension, 



will provide string corrections to the model. It would be interesting to estimate 
the effect of this using present local mass analysis. 

For further investigation it may also be interesting to analyze a black hole 
model proposed by Nojiri and Oda p0| . The model is described by SL(2,R)/U(1) 
gauged WZW model deformed by (1,1) operator and they claim that the black 
hole reaches a zero mass state and evaporates completely. 
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A. The Numerical Method 

The two dimensional surface (x + , x~) is represented by a two dimensional lattice 
with equal spacing in each direction. The lattice spacing in the x + direction is much 
narrower than that in the x~ direction. 

The equation ( |4.1|) can be symbolically represented by 

d_d + cf) = f{<p,p,d + <f)) , 

(Al) 

d-d + p = g(4>) d^d + <p . 

and for fixed x + the first equation can be seen as an ordinary differential equation 
of the first order for d + 4>. We use the fourth-order Runge-Kutta method to get d + 4> 
in the x~ direction. We have adopted the "Cubic Spline Interpolation" method to 
get values of an arbitrary x~ at fixed x + and d-<fi and d-p are computed with a 
fourth-order backward differential method. 

The values of <fi and p on the next x + = const, line are obtained from d + <p and 
d + p by linear extrapolation. 
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